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Partially invariant solution to (2 + 1)D shallow water equation is constructed and 
investigated. The solution describes an extension of a stripe, bounded by linear source 
and drain of fluid. Realizations of smooth flow and of hydraulic jump are possible. 
Particle trajectories and sonic characteristics on the obtained solution are calculated. 

1. Construction of the solution. Equations describing a thin layer of water flows 
over a flat bottom are observed. 

Ut + uu x + vu y + h x = 0, 

v t + uv x + Wy + hy = 0, (1.1) 
ht + (uh) x + (vh) y = 0. 

Here (u, v) is a velocity vector of a particle, h is a depth of fluid. The gravity acceleration 
is scaled to g = 1. Equations (II. lj) admit 9-dimensional Lie algebra of infinitesimal 
transformations L 9 [2], generated by operators (notations of [3] are adopted) 

X\ = d x , X2 = d y , X4 = td x + d u , X$ = td y + d v , 
X 9 = xd y - yd x + ud v - vd u , X w = d u 
X n = xd x + yd y + ud u + vd v + 2hd h , 

X12 = t 2 d t + txd x + tydy + (x — tu)d u + (y — tv)d v — Ithdh- 
X13 = 2td t + xd x + ydy - ud u - vd v -2hd h . 

We observe a partially invariant solution generated by Lie subalgebra 

N = {X U X 4 , X 10 + X 12 }CL 9 . 

Invariants of are the following functions 

A= ; I , V = vVt^Tl - 1\, H = h(t 2 + 1). 

According to the general algorithm [1] the representation of A r -partially invariant solution 
can be written as 

V = V(X), H = H(\), u = u{t, x, y). 

There are two invariant functions V and H, which are set to depend on the only invariant 
variable A, and one non-invariant function u, which depend on all independent variables. 
We substitute an expressions for invariants V and H in this representation, and then 
solve it with respect to unknown functions v and h. In what follows it is convenient to 
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represent function u as a function of (t, x, A) instead of (t, x, y). These give the following 
formulae: 

V(X)+tX , H(X) y 
u = u(t,x,X), v= \ J ^, h=— ^4, \= , . 1.2 

Substitution of the representation (11.21) into (II. ip gives 

u t + uu x + V{t 2 + l) _1 ltA = 0, 

l/y' + ^' = _A, (1.3) 

ra' + + (t 2 + ik - 1) = o. 

Here the lower indexes denote partial derivatives with respect to the corresponding vari- 
able; prime is the derivative with respect to A. The first and the third equations of system 
(11.3D form an overdetermined system II for non-invariant function u. 

At first, let us start from the subsystem II. Separation of variables in the third 
equation of fll.3p allows introduction of the new invariant function 

K(X) = -VH'/H - V. 

The non-invariant function u should satisfy the subsystem: 

V 

U t + UU X + ^-pyMA = 0, {t 2 + l) Ux ~t = K. (1.4) 

Integration of the second equation of ( 1 1.4ft gives 

u = x^±* + U(t,X). (1.5) 

Substitution of this representation into the first equation of (11.41) produces a compatibility 
condition in form of an equation for K: 

VK' + K 2 + 1 = 

and a relation for function U(t, A): 

{t 2 + l)U t + VU x + {K + t)U = 0. (1.6) 
Thus, the factor-system (\l.3\\ splits into the invariant subsystem 

VV + H' = -A, 

VK' + K 2 + 1 = 0, (1.7) 
VH' + HV = -KH. 

and equation (jl.6p for function U. These equations admit a discrete symmetry 

A -> -A, V -»■ -y, (1.8) 
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which follows from the admissible by equations (11 .11) transformation y — > — y, t> — > — 1>. 
Integration of the obtained nonlinear system of differential equations can be performed 
in terms of the new independent variable /z: 

dX Tr/ ,v f dX , 

rr v[x) ' " = Jvury (L9) 

With the new variable, the second equations of (II. 7p accurate to insufficient constant 
gives 

A" = -tan/i. (1.10) 
Taking ( 11.101) into account, one can integrate equation ( 11 .6D as 

/(/i — arctant) 
cos fiy/t 2 + 1 

(/ is an arbitrary function). Besides, system (11.71) has a first integral, which follows from 
its third equations: 

HV cos fi = m, m = const. (1-H) 

According to the discrete symmetry (11. 8ft it is enough to observe only positive values of 
V. Invariant depth H in non-negative by its physical definition, therefore, m cos/i > 0. 
As cos/i 7^ 0, the admissible set of values of functions and parameters can always be 
restricted to the following: 

AgR, ~\ < V i <\, V>0, A'(/i) > 0, m>0. (1.12) 

Finally, equations (11.71) possess a Bernoulli integral: 

V 2 + X 2 + 2H = b 2 (1.13) 

with an arbitrary constant b. This integral should be observed as an implicit (not resolved 
with respect to the derivative) ordinary differential equation for the dependence A(/i): 

(X') 2 + X 2 + -^— = b 2 . (1.14) 
v ' A' cos/i 

Thus, the solution is given by expressions ( 11.21) . (11.51) . in which functions V, H, K, and 
U can be determined from (jl.9p - (11.13p after the first-order ODE ( II . 14j) integration. This 
solution contain an arbitrary function / in the expression for the velocity component u. 
Equation ( II . 141) plays the key role in the further investigations. 

2. Properties of solutions of the key equation. For the analysis of equation 
( II. 141) it is convenient to observe a surface A, defined by the relation ( 11.141) in an extended 
space M 3 (/i, A, A') (see fig. [T]). Depending on the sign of the discriminant 

m 2 fb 2 -X 2 \ 3 , , 

cos z fi \ 3 / 

the cubic algebraic equation ( II . 14}) for the derivative A' possesses 1, 2 or 3 real roots. The 
surface A is correspondingly consists of three components A = Ai U A 2 U A 3 responsible 
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Figure 1: Representation of the implicit 
equation (11. 14ft smooth manifold 

A = Ai U A 2 U A 3 in an extended space 
M 3 (/x,A,A'). 



Figure 2: Integral curves of equation 
(11.141) on the surface A 2 UA 3 . The crimi- 
nant curve, separating A 2 and A 3 , divide 
integral curves into sub- and supercriti- 
cal parts. 



for each of three roots. The negative root X[ < exists for any value of parameters; the 
corresponding unlimited component of the manifold A is Ax (see fig. [T]). The remaining 
two components A 2 and A 3 corresponding to positive roots < A' 2 < A 3 are finite. They 
close up along the so-called criminant curve specified by an additional relation T> = 0. 
Physical meaning have only positive roots, since the negative root X[ = V < does not 
satisfy conditions (11.121) (the negative root corresponds to the negative water depth h). 
Surface A is woven from integral curves of equation (11.141) (see fig. In accordance to 
the general theory of equations, not resolved with respect to the highest derivative 
through each nonsingular point of the smooth surface A C IR 3 (yU, A,A') passes only one 
integral curve of the equation. The singular point could appear only at the criminant 
curve (= intersection of the cylinder T> = with the surface A 2 U A3). Solutions of the 
equation (11.141) in form A = A(/x) are obtained after the projection of three-dimensional 
integral curves into (/i, A)-plane. 

Thus, through each point of an area {(//, A) | T>{[i, A) < 0} pass two physically 
approved integral curves: one from the projection of A 2 , and another from the projection 
of A 3 . Inside of the area the integral curves have no singularities, therefore they can be 
found with any desired precision by the numerical integration of the equation (j!.14p after 
its explicit resolution for the derivative A'. The picture of the resulting integral curves 
is given in figure [3j In what follows the fluid flow is referred to as the subcritical one 
if the corresponding integral curve of equation fl!.14p lies entirely in A 2 surface; and as 
supercritical one if the integral curve lies in A3. Reasoning for this definition will be given 
below. 

Domain of each integral curve is the finite interval //*<//<//* (values and /x* 
are unique for each integral curve). At the boundaries of the domain (at the points of 
the discriminant curve T> — 0) subcritical and supercritical integral curves close up. At 
that, functions V(A), and H(X) are finite: 

V(K)= PM V3 ' V(X*) = ^W(X^=(^) 1/3 . (2.2) 

\COS/i*/ y COS /i J 

However, their derivatives go into infinity. Indeed, by the Chain Rule V'(X) = A"(//) / A'(/i). 
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Figure 3: Integral curves of equations ( 11. 14ft for m = b = 2. The domain of integral 
curves is limited by the discriminant curve. The solid (dashed) curves correspond to the 
larger (smaller) positive root A' of the cubic equation (11.141) . Bold integral curves satisfy 
boundary conditions //* = -1.450, A* = -2.341, /i* = 1.480, A* = -1.217 (dashed line), 
and /i* = -0.8, A* = -4.352, fi* = 1.251, A* = 3.833 (solid line). 




Figure 4: Dependencies ^(A) (solid line) and y/ H(X) (dashed line) for subcritical (left- 
hand figure) and supercritical (right-hand figure) flows. Calculations are made for the 
integral curves shown bold in figure [3J 



Equation ( 11. 14ft gives 

„ 2A' 2 + 2m sin/i(cos/i)~ 2 
= 3A' 2 + A 2 - 6 2 ' 

Positive roots of the cubic equation (11. 14j) coincide over the discriminant curve: A' 2 = A3. 
Hence, according to Viete's theorem, 3A' 2 2 = b 2 — A 2 , i.e. the denominator of the fraction 
is zero. Therefore, A"(//*) = A"(/x*) = 00. Graphics of dependencies V(A), and a/ H(X) 
for sub- and supercritical flows are given in figure HI 

3. Characteristics of the shallow water equations on the obtained solution. 
Let the characteristic surface of equations (11.11) be sought in the implicit form x(t, x, y) = 
const. There are three families of characteristics: contact ones 



Co ■ Xt + UXx + vx v = 
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Figure 5: Characteristic curves of the subcritical (left-hand figure) and supercritical 
(right-hand figure) flows, for integral curves of equation (j!.14p shown bold in figure El 



and two sonic ones 

C± : xt + uXx + vxy = ±VhyJxl + X 2 y - 

Let us investigate the set of characteristics, specified by the equation A = X{t). Sub- 
stitution of representations (|1.2p . (I1.5P gives equations for function \(t) along each of 
characteristics families in the form 



d\_V(X) dX _ V(X) ± y/H(X) 

In terms of new variables \i and r = arctant equations (13.11) can be integrated as 

A 

C-.T-^m. C ± :r-r =f m ^ 7m . (3.2) 
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Note, that for any chosen integral curve of equation (11.141) surfaces A = A* and A = A* in 
the space IR 3 (t, x, y) are sonic characteristics belonging to C_ family. They are enveloping 
surfaces for the rest of characteristics of C_ family, given by equation (13.21) . The typical 
characteristic curves (13.21) in (A, r)-plane are shown in figure [5l All the remaining ones are 
obtained by shifting curves in figure along Or-axis. In subcritical flow the disturbances 
velocity is higher then the fluid velocity, therefore, the disturbances are transmitted 
upstream reaching the left boundary A = A* of the domain of the solution. In the 
supercritical flow the situation is opposite: All the disturbances descend down the flow 
such that at the left boundary A = A* all characteristics are outgoing. 

4. Hydraulic jump. Let us show that the hydraulic jump may switch the super- 
critical flow to the subcritical one. Let the position of the hydraulic jump be fixed as 
A = Ao or, in physical variables, as y(t) = XoVt 2 + 1. Then the relative speed of the fluid 
motion is 

V{X )+tX X t V{X ) 
v n = v — y it) = . . = . 

w V^TT V^TT V^TT 
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As V > the state before the jump is to the left of the state after the jump. The 
conditions on the hydraulic jump (mass and momentum conservation laws) are [5]: 




(4.1) 



As usual, the square brackets denote the difference of the limiting values on the jump of 
the function in the states before and after the jump: [/] — fi — Substitution of the 
representation (11.21) allows rewriting of conditions (14.11) in the invariant form 



Let a = HV = mj cos fi. By virtue of the first condition of ( 14. 2ft we have [a] = 0. The 
second condition of ( 14.21) is convenient to rewrite as 



Function a/V 2 + 2V reaches its minimal value at V = <r 1//3 , i.e. at the boundary of the 
domain of the solution (see (12. 2p ). The flow is supercritical for V > a 1 / 3 , and subcritical 
for V < a 1 ' 3 . 

The jump stability condition is that in the state after the jump the fluid depth should 
be larger than the one in the state before the jump [5]. In terms of the solution under 
investigation this implies if 2 > Hi or, due to a conservation, V% < V\. This means that 
the state before the jump is supercritical, and the state after the jump is subcritical. This 
stability condition conforms with the Lax evolutionary condition for the strong disconti- 
nuity [6j. Indeed, graphics of characteristics in figure demonstrate, that incoming to the 
discontinuity A = Ao are 3 characteristics from the left-hand side and 1 characteristic from 
the right-hand side. The outgoing ones are two characteristics from the right-hand side 
of the jump. Thus, the number of outgoing characteristics is one less then the number of 
conditions (14. ip on the brake, which coincide with the Lax evolutionary condition. 

According to the characteristics equations (13.21) . function \l — t is Lagrange invariant, 
i.e. it conserves along particle's trajectory. As the trajectory of each particle is continuous 
at the hydraulic jump, values of \x are the same on both sides on the jump. By virtue 
of [a] = 0, constant m also conserves on the discontinuity: [m] = 0. Thus, the third 
condition of (14.21) is satisfied. The last condition of (14. 2p is equivalent to the demand of 
continuity on the jump of function / in expression for the velocity component U. 

Let the state "1" before the jump be known (i.e., constants m and bi are fixed and 
the solution is constructed according to formulae above), and the position of the jump 
A = Ao is specified. It is required to determine the state "2" after the jump, i.e. to 
select constant 6 2 7^ &i> such that the solution "2" conjugates with the solution "1" by 
conditions (14. 2p . 

The algorithm of solutions conjugation is the following. Note that the solution A(/i) 
of equation (11. 14ft is a monotonic function. Hence, value \x = \ii in the state "1" on 
the hydraulic jump can be uniquely found provided the jump position A = Ao is given. 
Known constants m and fii allow determination of a = mj cos/ii. Next, for known 
relative flow velocity V\ one can determine the velocity V2 after the jump from condition 
( 14. 3p . As a conserves on the brake, equation (11.141) gives constant b 2 in the state after the 



[HV] = 0, [2HV 2 + H 2 } = 0, [K] = 0, [U] = 0. 



(4.2) 



[a/V 2 + 2V] = 0. 



(4.3) 
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discontinuity. Let us show that obtained values of a, V2, and 62 satisfy condition ( 12 .11) of 
existence of physically meaningful solution of equation ( 11.141) . Indeed, 




Function V 2 + 2a /V reaches its minimal value 3<r 2//3 at V — cr 1 / 3 , which corresponds 
to the discriminant curve T> — 0. For all another values of V this function exceeds its 
minimal value, which guarantees the desired inequality T> < satisfaction. Finally, the 
solution behind the jump is constructed by finding the subcritical solution of equation 
( 11.14p with initial data A(/ii) = Ao- 

5. The general description of motion. Let us give an interpretation of the fluid 
motion governed by the obtained solution. For the fixed integral curve of equation fll . 14j) 
the motion takes place in a stripe, bounded by linear source and drain of fluid in positions 
y = \*Vt 2 + 1, and y = A*\A 2 + 1- Both lines are envelopes of the C_ characteristic 
family of equations ( 11. ip . Along Ox axis the stripe is not limited. However, due to the 
linear dependence of the velocity component u on the Cartesian coordinate x, the stripe 
should be limited with the use of, for example, two pistons. Thus, the flow takes place in 
a pool, whose two opposite walls are moving source and drain of fluid and another two 
walls are moving pistons. The free surface of fluid is convex upwards on subcritical flows, 
and is convex downwards on supercritical ones. The hydraulic jump along a moving line 
y = yoVt 2 + 1 is possible. In this case the convex downwards free surface of fluid switches 
stepwise on the jump to the convex upwards free surface. 

The particle trajectories are determined as the solution of Cauchy problem 

dx dy 

~dt =U) ~di = V, X ' i=0 = X °' y ' t=0 = Vo ' 
Substitution of the obtained solution and integration gives the following 

/i(A) — arctant — c\, x = (c 2 + /(ci) arctan t)\/t 2 + 1/ cos//. 

Here the dependence /i(A) is defined by the chosen integral curve of equation (11.141) . 
This dependence is monotonic, and, hence, is invertible at A* < A < A*. Let the value 
A = Ao = yo corresponds to /i = /io- Then equations of the particle trajectory, starting 
from the position (xo,yo) at t = are given by formulae 

x = (xo cos /io + /(/to) arctan t) Vt 2 + 1/ cos /i, /t(A) = /to + arctan t, y = \V t 2 + 1. 

(5.1) 

Lines y = const are frozen into the flow, i.e. they consist of the same particles at all 
moments of time. The fluid flow along Oy axis is fixed and is completely determined by 
the chosen integral curve of equation (11.141) . Fluid motion along Ox axes can be modified 
using the arbitrariness of function /. 

The completeness of the analytical investigation of the solution gives opportunity to 
utilize it for the testing numerical solvers of hyperbolic systems of equations with several 
independent variables [7]. 
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